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ABSTRACT 

Departures from complete frequency redistribution (CRD) in hydrogen lines are 
investigated for solar prominences. Partial redistribution effects (PRD) are found 
both in the wings (their already known lowering) and in the central part of the La 
line; a new feature is evidenced here: the partially coherent scattering in the 
near wings of the line leads to a double-peaked profile mirroring the incident 
solar radiation. With a low density model, we obtain a good agreement with 0S0 8 
observed profiles. On the contrary, the PRD computed LB profile (lower density, no 
reversal) departs from the observed one, a result which calls for more progress in 
terms of non l.t.e. transfer and modelling. 

1. INTRODUCTION 

During the last decade, several highly-sophisticated partial-redistribution techni- 
ques (PRD) have been applied to the analysis of the strong resonance lines formed 
in the solar chromosphere, plages or sunspots. Being low-density structures, solar 
prominences are also very good candidates for possible departures from complete 
redistribution (CRD), at least for most opaque lines like hydrogen La or LR.As 
reported by Vial (1982a), various high-resolution La and IB line profiles of a 
quiescent prominence have been recorded by the LPSP spectrometer on board of 0S0-8 
satellite. Applying the two-dimensional radiative transfer code of Mihalas et al. 
(1978), Vial (1982b) has simulated theoretical Let profiles (using a two-level atom 
and CRD) and made a comparison of their basic parameters with those corresponding 
to 0S0-8 data. The agreement was satisfactory, but no direct comparison between the 
shapes of the theoretical and measured profiles was presented in this study. It is 
the aim of the present paper to make possible a first comparison of 0S0-8 La and 
L B profiles with theoretical PRD computations in order to assess the real 
importance of departures from CRD. Moreover, we investigate the behaviour of the 
total CRD and PRD intensities of the hydrogen lines and compare our results with 
those of Milkey et al. (1979) and with the data presented by Heinzel and Vial 
(1983). 

2. PROMINENCE MODELS 

Since our principal interest is to estimate PRD effects on hydrogen lines, we start 
here with simple one-dimensional (ID) isothermal-isobaric prominence models. These 
models are of the (a - c) type of Heasley and Milkey (1976) or perhaps more 
realistic low-pressure ones of Heasley and Milkey (1978, 1983). The quiescent 
prominence as seen on the limb is represented by a vertically-standing slab of 
finite thickness with all model parameters being depth-independent. As basic input 
parameters we use the set (M, T, p, v fc , y) -see Tab. 1- where M is the total column 
mass along the line of sight, T is the kinetic temperature, p represents the total 
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gas pressure and v t characterizes the mean microturbulent velocity. The prominence 
plasma is assumed to be composed of hydrogen and helium atoms, the abundance ratio 
of helium to hydrogen being y = 0.1. 


MODEL 

M 0 

T 

p 

▼t 

y 


(g/cm 2 ) 

CK) 

(dyn/cm^) 

(km/sec) 


a3 

1 . 2-5 

7500 

0.065 

0 

0.1 

a5 

3 . 0-5 

7500 

0.065 

0 

0.1 

LP1 

1 . 2-5 

7500 

0.010 

0 

0.1 

LP2 

3 . 0-5 

7500 

0.020 

0 

0.1 

LP3 

3 . 0-5 

7500 

0.020 

6 

0.1 


Table 1. - Parameters of the prominence models used in this paper. 

Since the prominence slab is irradiated symmetrically from both sides by the 
incident solar radiation, we treat only one half of it. The corresponding second- 
order boundary conditions for the radiative transfer equation (RTE) are used in the 
same manner as in Mihalas et al. (1975), but are properly modified for the finite 
symmetrical slabs (see also Heasley and Mihalas, 1976). The surface boudary condi- 
tion requires the knowledge of the mean incident radiation field in all lines and 
continua. For Balmer lines we use the limb-darkened profiles of David (1961), Pa 
Pg and B a line profiles are taken from Zelenka (1976). All these profiles are 
subsequently adjusted to limb-darkened continuum level. Disc-center intensities of 
L a and Lg lines were obtained by 0S0-8 (see, e.g., Vial, 1982a). To obtain the 
intensity profiles of Ly and Lfi , we scaled the profile of LB by the ratios of 
total intensities Ly /Lg and Lfi /Lg , respectively. Integrated intensity of Lg 
is taken from OSO-8, Ly and L 5 from Vernazza and Reeves (1978). In order to obtain 
the most realistic mean-intensity profiles required in the surface boundary 
condition, we integrated numerically the above incident intensities, using the 
geometry and technique as discussed by Heinzel (1983). For all models we use the 
height 4000 km which corresponds to OSO-8 observations used here for a comparison. 
By the same procedure we also precomputed all fixed rates for optically thin lines. 
The continua are treated in two different ways: optically thick Lyman continuum 
transfer is solved in detail at several frequency points for which we specify the 
incident radiation temperatures equivalent to intensities measured on 0S0-6. For 
all subordinate continua we simply use one "mean" T ra j which defines the 
corresponding fixed rates (T^^ = 5480, 5900, 6200, 6400 K (without dilution) for 
Balmer, Paschen, Brackett and Pfund continua, respectively). 


3. FORMATION OF HYDROGEN LINES WITH PARTIAL REDISTRIBUTION 

The laboratory-frame redistribution function (LRF) used in the present computations 
takes generally the form (Heinzel and Hubeny, 1982) 

R (v», v) = y R i (v»,v ) + (1 -y ) R m (v, v ) (3.1) 

with R^ being either Rjj for the case of the resonance lines or Ry for the lines of 
subordinate series. While the functions Rjj and Ry contain certain degree of 
coherence of the scattering in the line wings, the function Rjh can be si®Ply 
replaced by complete redistribution in the observer’s frame. 
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The branching ratio y is defined as : 


y 


P J 

P j + Q E 


(3.2) 


where Pj is the total depopulation rate of the upper level j, while Qg represents 
the elastic collision rate. 

Strictly speaking, Eq. (3.1) has been derived assuming the impact approximation and 
isolated lines and, therefore, is not applicable to hydrogen lines. However, as 
demonstrated by Yelnik et al. (1981), the correct ARF (atomic-frame) for hydrogenic 
La line is formally the same as that corresponding to LRF (3.1), but the elastic 
collision rate Qg is now frequency-dependent. Consequently, the atomic profile 
function is no longer Lorentzian in the wings and, moreover, the branching ratio 
is also frequency-dependent (Qg in (3.1) is replaced by Qg ( v )). 

For LB we apply the same procedure with Qg (v) taken from Gouttebroze et al. 
(1978). For higher members of the Lyman series, we simply use CRD. Moreover, 
Heinzel (1983) has shown that for quiescent prominences observed on the limb, the 
subordinate lines exhibit negligible departures from CRD when applying LRF R y . 
Anticipating this result, we use in the present computations CRD instead of 
complicated LRF (3.1) for all subordinate lines. 

As a basic procedure to solve the non-LTE transfer problem with PRD we use a 
combination of an equivalent-two-level-atom approach (ETA) and partial-lineari- 
zation technique. In the frame of the so-called standard PRD problem (see below), 
ETA was modified by Hubeny ( 1 985) to account directly for PRD 
in one line at a time. The resulting source function is 


S (v) =- 


(v) -7) + n 


1 + 


(3.3) 


where J is the integrated mean intensity and Rjj (v) is the redistribution integral. 
The parameters n and ev in (3.3) have the similar meaning as in standard ETA (for 
details see Hubeny, 1985). A very important parameter appearing in Eq. (3.3) is the 
parameter A 



(3.4) 


where Aj^ is the spontaneous emission rate for the given transition. For strong 
resonance lines like La or LB , A in form (3.4) reflects the coherence properties 
of the transition i 4 — > j. Equations (3.3) and (3*4) are valid if all multilevel 
cross-redistributions lead to a natural population of the upper level j - this is 
exactly what we call the standard PRD problem, i.e. we consider here no PRD 
multilevel interlocking. 

Simply speaking, y alone is the branching ratio derived for a two-level-atom, while 
A accounts for multilevel processes. From (3.4) we get for L« A = y since P 2 = 
A 21 . On the other hand, for LB we have P 3 = A 31 + A 30 30 that A = 0.56y . Since the 
ratioyfor both La and LB is nearly equal to unity (for low-density media), the 
importance of the redistribution term R-q (v) in (3*3) is significantly reduced for 
L B due to the scattering process 2 -> 3 -> 1 which is assumed to contribute to the 
emission in L {5 as a CRD process. 

The numerical procedure for solving the full non-LTE problem is the following: the 
basic iteration loop uses the linearization scheme as described by Mihalas et al. 
(1975) to solve simultaneously RTE (by Feautrier method), equations of statistical 


157 


equilibrium, particle and charge conservation equations. Hydrostatic equilibrium is 
treated iteratively, the temperature structure is given from the model. We use 5- 
level hydrogen atom with continuum and explicitely treat all Lyman lines, Lyman 
continuum and Ha line. Other transitions are assumed to be optically thin and their 
radiative rates are fixed by the external radiation field. The linearization loop 
ensures the global convergence of this highly non-linear problem. Between each two 
linearization iterations , we perform several ETA iterations , accounting directly 
for PRD in La and LB lines. In this way, we obtain fast convergence (within 2-3 
linearization iterations) in all 50 depths and for all transitions and level 
populations. We use most important opacity sources for hydrogen spectrum formation 
in low-density prominences, both ETA and linearization schemes use variable 
Eddington factors (Auer and Mihalas, 1970) to minimize the computer time. Inelastic 
collisional rates for hydrogen are taken from Mihalas et al. (1975). The redistri- 
bution matrices are used as depth-independent since Pj>> Qe anc * the electron and 
proton densities vary only slowly with depth. On the other hand, the branching 
ratio y(v) is computed for all depths. 

4. DISCUSSION OF NUMERICAL RESULTS AND COMPARISON WITH OBSERVATIONS 

In Fig. 1 we display typical La profiles for both PRD and CRD approaches (for the 
model LP1). 

Drastic differences between CRD and PRD profiles can be explained in the following 
manner. Around line center, PRD behaves similarly as CRD and the emergent radiation 
is saturated at the level of diluted incident chromospheric intensity. However, 
when the scattering takes place in the near wings, we can observe an interesting 
effect: 



Fig. 1 - La intensities emergent from a prominence slab (model LP1). Both PRD (-) 

and CRD ( ) profiles are displayed for three outgoing angles with y = 1., 

0.6 and 0.2 (broader profiles correspond to lower values ofy ). 
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for CRD, the photons absorbed in this region 
can freely diffuse into the wings which 
become broad and intense; PRD is, on the 
other hand, represented by partially- 
coherent scattering in the wings so that the 
incident chromospheric L peaks (centered at 
AA=i0.2 R) are partly "reproduced" for the 
PRD case. The photons absorbed in the core 
or in the near wings cannot freely penetrate 
into the wings as in the CRD case and, as a 
consequence, far-wing intensities are suffi- 
ciently below the CRD values. Very low PRD 
wings follow from the well-known behaviour 
of the redistribution function Rjj* 

Partial reproduction of the strong incident 
La peaks has the following consequences. 
First, due to quasi-coherent penetration of 
these peak photons into the slab center, we 
arrive at higher excitation of hydrogen, 
i.e. the populations of the excited levels 
are higher for PRD. As a consequence, the 
lines of subordinate series are more intense 
for PRD than for CRD, which is demonstrated 
for Ha line in Fig. 2. 

In this figure, La integrated intensities 
are nearly the same for both PRD and CRD 
(for PRD, the low intensities in the far 
wings are compensated by the two peaks), but 
Ha is shifted in PRD case. In this way we 
can also explain the results obtained by 


Fig. 2 - Dependence of the ratio E(Ha )/E(La ) on E(Ha ), E being the integrated 
intensity in units ergs s" 1 cm -2 sr" 1 . Hatched region contains the observed ratios 
taken from Heinzel and Vial (1983), dashed lines correspond to CRD intensity-ratios 
as obtained by Heasley and Milkey (1976) for the models indicated in parentheses 
(labeled by kinetic temperature and total gas pressure). Our theoretical ratios for 
models a3 and a5 are to be compared with squares on the dashed curve (7500, 
0.065). Our CRD (o) and PRD (•) computations are presented for P = 1. 

Milkey et al. (1979): their CRD intensity ratios are shown in Fig. 2 as two dashed 
lines, al - a5 models lie on the curve (7500, 0.065). For models a4 - a5, this 
curve is significantly below the (hatched) region of observed values as taken from 
Heinzel and Vial (1983) - this discrepancy is caused probably by adopting unrea- 
listic incident La radiation field. The theoretical L a profiles of Milkey et al. 
(1979) exhibit only the second PRD-feature, i.e. the low wing intensities. These 
authors obtained no difference between CRD and PRD for the line core and near wings 
and, subsequently, their integrated La intensities are a factor two lower for PRD 
as compared with CRD which leads to an apparent agreement with Skylab data. Further 
more, we note that LP1 is the model with lower M and, therefore, for »»= 1. and 0.6 
the incident peak is not fully reproduced due to lower optical thickness of L a (see 
Fig. 1). However, for thicker models like LP2 (Fig. 3), LP3 or a5, the incident 
peak is nioely reproduced and its position agrees well with that for the original 
peak. It seems that the observed L« peak-to— peak distance (PTPD) in quiescent 
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prominences is simply a measure of partially-coherent reproduction of the incident 
chromospheric peaks and the actual value of PTPD depends on the coherence ratio! 
and on the total opacity in La . 

Finally, PRD leads also to higher electron densities inside the slab, simply due to 
higher excitation of the second level which controls the photoionization by the 
Balmer continuum. 

For a direct comparison with 0S0-8 observations (Vial, 1982a) we used the model 
LP2. 



- 0.99 - 0.15 - 0.50 - 0.25 0 . 0.25 0.50 0.15 0.99 


WAVELENGTH ( A ) 


Fig. 3 - Comparison of the observed Lx profile (full line) with theoretical PRD and 
CRD intensities for model LP2 (y = 0.6). For details and discussion see 
the text. 

In Fig. 3 we display the observed La profile, together with corresponding theore- 
tical PRD and CRD profiles. PRD profile has been smoothed to account schematically 
for the finite spectrometer resolution (0.02 A) and for small line-of-sight velo- 
cities of the order of 5 km/s (see Jensen, 1982). While CRD profile is in no case 
capable of explaining the observed shape, our PRD La profile leads to a reasonable 
fit. However, there are certain differencies to be discussed. First, the difference 
in the central dip can be due to either small inaccuracies in the geocoronal- 
absorption removing, or simply to inappropriate value of the incident radiation 
used in the computations. Secondly, the observed peaks are somewhat higher than 
computed. In fact, for thicker models like LP2 we arrived at rather flat peaks and 
these can be slightly modified when the deconvolution procedure is applied (this is 
also valid for the flat core). Finally, a wavelike character of the far wings is 
due to deconvolution so that we can rely on the intensities up to approximately + 
0.5 A from the line center. As we have found from the computations, the central 
intensity of Lx is determined mainly by the value of the diluted incident radiation 
and is insensitive to variations of the model parameters. On the other hand, the 
peak value is somewhat lowered for non-zero v fc (model LP3). 

The profile parameters, corresponding to Fig. 3, are summarized in Table 2. We see 
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that our PRD values reproduce quite well the observations, possibly except for the 
total intensity corrected for the filling factor (note that the actual observed 
value of E can be lower with respect to data-reduction effects discussed above). 
The 2D total intensity is too low even when CRD is used. Note also a satisfactory 
agreement between H o intensities. Finally, the electron density as obtained for LP2 
model (n e = 7 x 10* cm“^) is in good agreement with the values recently obtained 
from polarimetric observations (Bommier, present volume). The geometrical thickness 
of the prominence is somewhat large - 11 500 km -, but this corresponds to a low 
gas pressure 0.02 dyn/cm . 



La 

LB 

Ha 


LP2 

LP2 

2D 

OBS. 

LP2 

LP2 

OBS. 

LP2 

LP2 OBS. 

PARAMETERS 

PRD 

CRD 

CRD 

FPO 

PRD 

CRD 

FPO 

PRD 

CRD 

FWHM (A) 

0.76 



0.74 

0.26 

0.23 

0.61 



PTPD (A) 

0.37 

0.23 

0.32 

0.37 

(no reversal) 

0.33 



I(peak)/I Q 

1.25 

1.06 

1.3 

1.33 

(no reversal) 

1.85 



Ei (xIO 4 ) 

2.70 

3.03 

1.61 

2.84 

0.007 

0.005 

0.044 

7.51 

6.76 5—9. 

E 2 (x10 4 ) 




3.55 



0.055 



PTPD - peak-to-peak distance - I Q - central intensity (prominence) 

Ei (ergs s -1 cm -2 sr -1 ) - integrated intensity 

E 2 - as E-| but corrected for filling factor (see Vial, 1982a) 


Table 2. 


Comparison of basic profile parameters of L« , Lg and H a » LP2 is the model from 
Table 1, 2D model corresponds to the computations of Vial (1982b), FP0 means "first 
part of orbit" (see Vial, 1982a). Ha intensity was estimated using 0S0-8 Call 
intensities and the relation between Call and H a as given by Stellmacher (1979). La 
parameters correspond to the profiles in Fig. 3. Theoretical Lq profiles for the 
model LP2 have been corrected for the actual value of the calibration disc- center 


intensity (6.5 x 


10 4 ergs s” 1 


cm 


-2 „„-1 


sr -1 ), 


for LP2 we used the values for y = 0.6. 


As we have found during the course of our investigation, Lp line presents, on the 
other hand, a complicated and so far unresolved problem. Looking at Table 2, we see 
immediately that the basic Lg profile-parameters differ significantly from the 
observed ones and, moreover, the theoretical integrated intensity is also a factor 
5-10 lower than the observed value (for both CRD and PRD). This surprising result 
needs further verifications and more sophisticated calculations are to be done to 
resolve this discrepancy. 

Our knowledge of the L g radiation field inside the prominence is also essential for 
estimating the degree of Ha line polarization (from which we can deduce the 
magnetic-field topology via the Hanle effect -see Bommier, present volume-) and, 
therefore, PRD-interlocking as well as various multidimensional effects are to be 
properly accounted for in the future computations. 
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5. CONCLUSIONS 


We have demonstrated the most important effects of PRD in the formation of hydrogen 
La line in quiescent prominences. Except for the well-known lowering of the far 
wings as compared to CRD, we have found a new PRD feature for L<, prominence line: 
strong peaks of the incident solar radiation are partly reproduced due to quasi- 
coherent scattering in the near wings of L« and, as a consequence, we obtain typical 
prominence La profiles with significant central reversal and PTPD comparable to the 
solar one. All features we have found are actually observed as our comparison with 
0S0-8 data shows. However, for Lg line we arrived at a substantial discrepancy 
between the theoretical and observed intensities. Having estimated the differential 
PRD effects for La , future computations should try to establish the proper 
interplay between PRD and multidimensional radiation-transport influences on the 
prominence spectral diagnostics. Our present report should serve as a starting point 
in developing more sophisticated non-LTE techniques which will be capable of 
explaining new data from UVSP/SMM or other planned space experiments. 
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